CERN-TH/2001-229 
August, 2001 



Dressed gluon exponentiation^ 



Einan Gardi 

TH Division, CERN, CH-1211 Geneva 23, Switzerland 



Abstract: Perturbative and non-perturbative aspects of differential cross-sections close 
to a kinematic threshold are studied applying "dressed gluon exponentiation" (DGE). 
The factorization property of soft and coUinear gluon radiation is demonstrated using 
the light-cone axial gauge: it is shown that the singular part of the squared matrix 
element for the emission of an off-shell gluon off a nearly on-shell quark is universal. We 
derive a generalized splitting function that describes the emission probability and show 
how Sudakov logs emerge from the phase-space boundary where the gluon transverse 
momentum vanishes. Both soft and coUinear logs associated with a single dressed gluon 
are computed through a single integral over the running-coupling to any logarithmic 
accuracy. The result then serves as the kernel for exponentiation. The divergence of the 
perturbative series in the exponent indicates specific non-perturbative corrections. We 
identify two classes of observables according to whether the radiation is from an initial- 
state quark, as in the Drell-Yan process, or a final-state quark, forming a jet with a 
constrained invariant mass, as in fragmentation functions, event-shape variables and deep 
inelastic structure functions. 
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1 Introduction 



Many observables in QCD receive large corrections from soft and collinear gluon radiation. 
Contrary to virtual corrections, real gluon emission has a restricted phase space that 
depends on the kinematics. In particular, close to a kinematic threshold, a differential 
cross-section is dictated by soft and collinear gluon radiation. Classical examples are 
event-shape variables close to the two-jet limit, and fragmentation functions, deep inelastic 
scattering structure functions and Drell-Yan production for x — > 1. 

Differential cross-sections close to threshold have been studied extensively in pertur- 
bation theory, where the phenomenon is reflected in Sudakov logs L = In 1/(1 — x). The 
presence of large logs in the perturbative coefficients limits the applicability of fixed-order 
calculations. It is a universally acknowledged conjecture that perturbation theory can be 
improved by resumming the leading logs to all orders in the coupling ||l|-[^. 

The resummation of Sudakov logs is based on the general property of factorization. 
Factorization of soft photons (gluons) goes back to the Low theorem P, |^, which states 
that the leading, 0{l/w) {w is the photon energy), and next-to- leading, 0{w^), terms in 
the radiative amplitude in a generic process can be expressed in terms of the radiationless 
amplitude. As a consequence, in the soft approximation the amplitude can be written 
as a product of a universal bremsstrahlung factor times the amplitude that describes the 
hard process. Under certain conditions, the bremsstrahlung factorization theorem can be 
extended to include hard collinear radiation with small transverse momentum kt (see 
also i-i and gM)- 

An important consequence of factorization is that the log-enhanced part of the per- 
turbative expansion of differential cross-sections exponentiates. This means that a leading 
order calculation of the single parton branching probability is sufficient to generate the 
leading contributions to any order in perturbation theory. 



In single-scale observables renormalization group and physical considerations like [|lT 
allow one to choose the renormalization scale /i^ such that unnecessarily large logs are 
avoided. In differential cross-sections there are several physical scales, e.g. the hard 
scale and Q'^{1 — x), so large logs are inherent. The choice of the renormalization 
point as the hardest scale /i^ ~ will make the coupling small. However, for x — 
1 the cross-section is determined by multiple soft and collinear emission with tj^ical 
momenta that are much smaller than Q. Using OsiQ"^) as an expansion parameter will 
force the coefficients to be large, reflecting the presence of the lower scale Q^(l — x). As a 
consequence the resummation of the differential cross-section is strictly restricted to the 
range Q^(l — x) > A^ where A is the scale at which the running coupling diverges. In 
addition, within the perturbative domain, perturbative and non-perturbative corrections 
are large. 

For fragmentation functions and deep inelastic structure functions, which are both 
governed by collinear emission, the perturbative domain is Q^(l — x) > A^. In event-shape 
distribution, and in Drell-Yan production close to threshold, large angle soft emission 
enforces a more stringent restriction, e.g. Q{1 — T) > A in the case of the thrust (T). In 
terms of the coupling, the largest possible domain of applicability of Sudakov resummation 
is therefore AL < .^max, where A = /5oas(<5^)/7r and ^max = 1 or 1/2 in the two cases 
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mentioned. 

It is well known that perturbative coefficients in QCD increase factorially at large 
orders. The dominant source of divergence are diagrams that are related to the running 
coupling, the so-called "renormalons" . Due to renormalons, the applicability of fixed- 
order perturbation theory may be more restricted than a priori expected. The problem 
can be addressed in perturbation theory by the resummation of the running coupling 
effects, using the analogy with the Abelian theory • 



It is important to realize that a perturbative treatment of running coupling effects 
is not sufficient: infrared renormalons make the resummation procedure ambiguous at 
power accuracy. For some inclusive observables, such as the total cross-section in e^e~ 
annihilation, power corrections are small ~ This can be established by an operator 

product expansion or by the analysis of infrared renormalon ambiguity |15|, [l^, 0. On 
the other hand, in differential cross-sections, the leading corrections are determined by 
the lower scale, i.e. ~ 1/(Q^(1 — x))" (or, for event-shape variables such as the thrust, 
1/{Q{1 — T)Y). Consequently also the power-correction expansion is expected to break 
down for Q'^{1 — x) ^ (or for Q(l — T) ^ A in the thrust case). The positive side is 
that the factorization property holds beyond the perturbative level. It is natural to expect 
that like the perturbative series also the power-correction series exponentiates close to the 
kinematic threshold [|T7| -p3|. 

The impact of renormalons on differential cross-sections is both interesting in principle 



and important in practice. The analysis of the thrust distribution e^e annihilation in |]23 



showed that due to renormalons, sub-leading logs are factorially enhanced with respect to 
the leading logs. This observation is completely general: this is the way in which sensitiv- 
ity to lower scales is reflected in the perturbative coefficients. The immediate consequence 
is that resummation based on a fixed logarithmic accuracy has a limited range of validity, 
LA <^ ^max- A quantitative treatment in a wider range (still LA < ^max) requires to sum 
all the factorially enhanced sub-leading logs and include the associated power corrections 
in the exponent. This is achieved by Dressed Gluon Exponentiation (DGE) [^]. 

DGE is designed to deal with both Sudakov logs and renormalons. It is based on 
exponentiating the entire perturbative series (of log-enhanced terms) that describes the 
single gluon emission close to the threshold. At a difference from the standard parton 
cascade approach (see e.g. 0, 0), the single gluon emission probability is calculated to 
any order in perturbation theory (in the large-/3o limit) and to any logarithmic accuracy, 
thus incorporating the factorial enhancement of sub-leading logs. In addition to the 
improved perturbative treatment, by identifying the divergence of the exponent one can 
introduce parametrization of the relevant non-perturbative corrections. 

A good example where renormalon analysis is crucial in identifying the leading power 
corrections is the case of Drell-Yan production |T^, p5[-|2^. Here, kinematic consider- 



ations suggest that as in the thrust distribution, soft gluon emission at large angles will 
result in a 1/(5(1 — x) correction. Renormalon analysis ||26|, ^ shows that the leading 
correction is in fact ~ 1/Q'^{1 — x)"^. This difference between the case of event-shape 
variables and that of Drell-Yan is very intriguing and we shall return to discuss it below. 
In both and [Q, the starting point for the perturbative calculation was that of 



the dispersive approach in renormalon calculus 03, E^] (see also 113 )• evaluation of 
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the single dressed gluon (SDG) cross-section based on the exact matrix element and the 
exact off-shell gluon phase space. On the other hand, the factorization property and 
the resulting Sudakov log resummation methodology (see e.g. @, |, ^, |2^) are based on 
keeping the gluon on-shell and making specific kinematic approximations that single out 
the singular parts of the gluon emission matrix element. Although integrals over the 
running coupling have a central role in both cases, the scale of the coupling in the former 
is the gluon virtuality, whereas in the latter it is its transverse momentum. Thus at first 
sight, the two methodologies seem mutually exclusive. In this paper we show how they 
can be bridged over. We identify the appropriate approximation required to single out 
the singular part of the matrix element in the case of an off-shell gluon emission. Next, we 
show that the factorization property applies to off-shell gluons in the same way it applies 
to on-shell gluons. This allows us to fully exploit the power of factorization: the entire 
renormalon sum exponentiates. 

The calculation of the exponent is, of course, approximate. However, instead of work- 
ing at a fixed logarithmic accuracy, our approximation is based on the Abelian large-/3o 
limit. The usual reason to consider this limit is that it identifies the dependence on the 
scale of the running coupling and thus also the infrared sensitivity. In the context of 
exponentiation there is an additional motivation: the large-/3o limit distinguishes single 
gluon emission from multiple emission. Since the exponentiation kernel is primarily as- 
sociated with a single gluon emission, it is natural to evaluate it using this limit. The 
first step beyond the large-/?o limit is taken in an unambiguous way by identifying the 
running coupling with the "gluon bremsstrahlung" effective charge [^,2^. Using this 
coupling and the two-loop renormalization group equation is sufficient to guarantee that 
the DGE resums all the logs (i.e. not only the large-/?o logs) up to next-to-leading log- 
arithmic accuracy (NLL). Contrary to a fixed logarithmic accuracy resummation, DGE 
is free of renormalization-scale dependence. Clearly, some residual scale dependence ap- 
pears beyond the logarithmic approximation upon matching the exponentiated result with 
fixed-order calculations [^. However, in the threshold region, when the logarithmically 
enhanced terms dominate, this effect is negligible. 

The logarithmically enhanced cross-section can be calculated from the singular terms 
in the squared matrix element. Using the light-cone axial gauge we demonstrate that as 
far as the singular terms are concerned, the emission probability is process-independent. 
The squared matrix element factorizes into a hard part times a function that depends only 
on the gluon momentum. The latter is interpreted as a generalized splitting function. The 
factorization of the cross-section depends on the phase space. Here we consider the case 
where all energetic particles move in two light-cone directions. Such light-cone kinematics 
arises in different cases: in e"'"e~ these are the two opposite directions of the quark jets, 
in Drell-Yan production near x = 1 these are the directions of the two incoming quarks, 
and in deep inelastic scattering near xb] = 1 one light-cone direction is the direction of 
the incoming quark, the other that of the outgoing quark. In each of these cases we show 
how the singularity of the matrix element translates into Sudakov logs through the phase- 
space integration. It is crucial that the logs emerge only from the limit in which the gluon 
transverse momentum vanishes (it can be soft or collinear). As a consequence, the relevant 
terms in the characteristic function can be easily computed. The resulting characteristic 
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function is integrated with the running couphng yielding a resummed cross-section that 
contains all the logarithmic terms associated with a single dressed gluon. This is the 
kernel for DGE. 

The outline of this paper is as follows. In section 2 we use the light-cone gauge to 
analyse the generic scenario of an off-shell gluon bremsstrahlung in light-cone kinematics. 
We identify the approximation needed to isolate the singular parts of the squared matrix 
element and show that in this approximation the radiative squared matrix element can be 
written in a factorized form in which a universal bremsstrahlung factor multiplies the non- 
radiative squared matrix clement. In section 3 we turn to discuss factorization at the level 
of the cross-section. Considering first the case of final-state radiation in e+e" annihilation 
we derive a generalized splitting function formula for the case of an off-shell gluon. In 
section 4 we illustrate the advantage of the axial gauge in this context. Using the case of 
the quark fragmentation function as an example, sections 5 and 6 describe the two steps 
in promoting the leading order calculation to a resummed one: in section 5 we show how 
phase-space integration can be performed within the approximation considered to obtain 
the relevant part of the SDG characteristic function and then use the dispersive approach 
to dress the gluon through an integral over the running coupling. In section 6 we show how 
multiple emission can be taken into account by exponentiation of the SDG cross-section. 
Sections 7 though 9 summarize the application of the method to other observables: jet 
mass (or thrust) in e~^e~, coefficient functions in deep inelastic scattering and Drell- 
Yan production. Common features and differences between the different examples are 
emphasized. Section 10 contains some concluding remarks. 



2 Factorization of the squared matrix element 

Sudakov logs appear because of soft or coUinear gluon radiation off nearly on-shell partons. 

The emitting parton can be a gluon or a quark. In the standard treatment of parton 
cascades, each possible branching is assigned a probability, which is a function of the 
longitudinal momentum fraction. Here we develop a similar probabilistic formalism for 
the case of an off-shell gluon emission off a quark. Gluon splitting will be taken into 
account in a different way, through an integral over the running coupling. 

In a generic QCD process, consider gluon emission off an outgoing quark, assuming 
that the quark is on-shell after the emission = rrig ~ 0, while the gluon possesses a 
time-like virtuality k"^ = im?. The amplitude for the emission off this particular quark is 

^-9st'' u^'^iph^ijf + l^)Mr, (1) 

where is a colour matrix, e^* is the gluon polarization vector and M.r represents the 
rest of the process. The enhancement of soft or coUinear emission is associated with the 
singularity of the quark propagator 1 / (k which is only partially compensated by the 

terms in the numerator. Factorization of the cross-section is based on the fact that the 
singularity emerges just from this diagram. The squared matrix element corresponding 
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Figure 1: Gluon bremsstrahlung off an outgoing quark 



to (|T|) summed over the gluon and quark spins is 

E MM^ = fE^r^;) Tr {(2/ + + l^)MrM;} (2) 



where Y.s'^^^Kp)'^^^Kp) was replaced by ^. Since the gluon attaches also to other parti- 
cles with an amplitude A4, the full calculation of the cross-section should be based on 
{Ad + + Aiy rather than on (|^). However, the explicit terms in @ capture the 

entire singularity if interference terms of the form M.Ai'^ are non-singular. This depends 
on the gauge choice X^a ^t* = d^iu- 

Defining two light-like vectors hj p = 0, 0) and p = (0,p_, 0), where = = 0, 
we choose the light-cone axial gauge where = 0, which is a physical gauge: g^'^d^u = —2. 
The gluon propagator is given by 

df,u = -g^.u + • (3) 

Let us consider now light-cone kinematics. Particles moving in the "— " direction couple 
to the gauge field through J_y4+, where J_ is an effective current. Since in this gauge 
Aj^ vanishes, these particles do not interact with the gauge field. Therefore, if we assume 
that the other coloured particles move close enough to the "— " direction, it is sufficient to 
consider the amplitude in (|1|). Put differently, in this gauge (0) captures the entire singu- 
larity and the interference terms M.M.'^ are non-singular. This point will be demonstrated 
explicitly below. 

Note that the assumption we made here is rather strong: the invariant mass of the 
particles moving in the "— " direction must he negligible. This assumption is relevant in 
many applications including event-shape variables close to the two-jet limit, and fragmen- 
tation functions, deep inelastic scattering structure functions and Drell-Yan production 
for X — i> 1. On the other hand, this assumption is inadequate for a generic jet observable 
in hadron-hadron collisions where both initial- and final-state radiation are important, 
nor for the case of three jets in e"'"e~ annihilation. 
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Let us now calculate the squared matrix element (0) in the light-cone gauge. The 
traces corresponding to the two terms in the propagator are 



ri = 2 Tr 



r2 = 



- lj:)MrMr^] 
^ [Tr + l^)m{T^ + ^)MrMr^} + Tr 



(4) 



kp 



{4pp + 2kp) Tr {^MrMr^} + 2ppTr {^MrMr'^} - 2pkTT [^MrMr^} 



Collecting the terms r = ri + r2 one obtains 

jk + pf 

kp 
{k + pf 



-2m' 



Apk + 2pp 



kp 



{App + 2kp) 



-2pk{k + pf 
kp 



(5) 



Introducing the Sudakov decomposition of A; = kt) = Pp+ap+kt, and defining 

the following dimensionless parameters 



P = k^/p^ = 2kp/2pp 

a = k^/p^ = 2kp/2pp 

7 = kl/2pp 

A = /2pp = (2k+k- — k^)/2pp = a(3 — 7 



(6) 



the trace becomes 

r = 2ppTi[^MrMr^] 

+ 2ppTT [l/^MrMr^] 



-2A + 2^4^(2 + 



2a + 2 



a + \ 



+ 2ppTT{:^MrMr^} 



-2a 



a + \ 



(7) 



Sudakov logs in the differential cross-section are associated with the region where 
the denominator of the quark propagator {p + k)"^ = 2pk + fc^ = 2pp {a + A) is small. 
Therefore, from now on we assume a + A -C 1. Since both a and A are positive, both 
must be small. The gluon on-shell condition ajS = 7 + A then implies that the transverse 
momentum fraction 7 is small, at least with respect to (3. We are therefore interested in 
the contribution from the region where a, A and 7 are all small, but no specific hierarchy 
between them is imposed. This approximation is adequate to treat both the soft and the 
collinear regions: f3 is small in the soft limit, but it is of order 1 or larger in the collinear 
(non-soft) limit. 

In this approximation one can simplify (|^) by 

• replacing Tr |^7VlrA^r^} by /STr |]/A^r-^r^}- The other components of k have 
additional suppression by one of the small parameters. 
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ignoring the terms proportional to Tr |^7Vlj.-^r^|- Taking p_ of the order of p+, 
these are exphcitly suppressed by a or A, with respect to the relevant terms. 



Then the following expression for the trace r is obtained 

(3^ + 2p + 2 



-X{l + /3) + {a + X) 



/3 



(8) 



and, finally, the squared matrix element (|) is 

-A ,^ 1 (3^ + 2(3 + 2 



J2MM^ = 2CFgl 



a + X 



P 



2pp 



(9) 



This demonstrates that as far as the a + X — > singular terms are concerned, the 
squared matrix element is factorized into a process-independent factor, corresponding to 
the emission of an off-shell gluon, times Tr |]/A^,,A^r^}; which is the expression for the 
squared matrix element in case of no gluon emission. The latter is, of course, process- 
dependent. 

In the case of gluon radiation off an incoming quark, similar considerations yield 



1 



J2MM^ = 2CFgl 



-X 



{a-xy 



1 /3^-2(3 + 2 



a — X 



2pp 



(10) 



where the factor 1/2 accounts for averaging over the quark spin. This expression can be 
obtained directly from (P) by taking a — > —a and (3 — > —j3, corresponding to inverting 
the momenta p and p. Note that in ([T0|), as in (^, the Sudakov parameters a, (3, 7 and 
A are all positive. 



3 Gluon emission probability 

To show factorization of the physical cross-section, one has to consider the phase space. 
Suppose that the external particle momenta are denoted by p^ and p^ for incoming and 
outgoing particles, respectively. The cross-section for p\ + — ^ J2n pi + k is 

1 _ _ r d'^p^ 1 

4EaEb\va-vb\ ny J2^^ ^^^^ 

Because of momentum conservation, the gluon is correlated with the other final state 
particles. Nevertheless, kinematic factorization can be established for a simple enough 
phase space in the appropriate Mellin or Laplace space. The details of factorization and 
the resulting exponentiation formula depend on the observable (see e.g. 0, H, ^). 

As an example we shall consider here (and in sections 4 through 7) the case of e~^e~ 
annihilation into hadrons. It is convenient to identify the vectors that define the and 
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"— " directions p and p as the final momenta of tlie two primary quarks. Working in tlie 
approximation of independent emission, we can calculate the gluon emission probability 
as if there are only three outgoing particles p, p and k. Multiple emission results, in this 
approximation, in straightforward exponentiation, as described in section 6. The centre- 
of-mass energy squared is g^, where q = p\+ p^^. Note that it is not given by 2pp: the 
contribution of gluon momentum cannot be neglected in the collinear limit. 

Assuming that factorization occurs, we proceed by performing the integration over p 
and p in the trivial case where Tr |]/A^r-^r^} is proportional to 2pp. The three-particle 
phase-space integral is then 

^ ^ dadl3- — — ~ / dadp . (12) 



1287r3 7 ^(l + a + /5 + A)3 1287r3 7 (1 + /?) 

while the two particle phase-space integral is 

Starting with the squared matrix element (^ integrated over the three-particle phase 
space, we now identify the gluon emission probability as the factor multiplying the cross- 
section for no gluon emission (which is Tr |]/A^f.A^r^| integrated over the two-particle 
phase space). Changing variables from A to 

^ A {2pp) ^ \ ^ 

g2 i + a + p + X l + ^ ' 

guarantees that the integration is performed for fixed and Q^. The gluon emission 
probability is given by 



dP- ^^"^ 



27r 



-A 1 P^^20^2 



{a + Xy a + A (3 



dad(3dX6(X-e(l + p)) , , 

■ 



Equation ( |T5| ) can be interpretated as a generalized splitting function for the off-shell gluon 
case. The standard splitting function can be obtained in the limit A = 0: substituting 
z = 1/(1-1-/3) as the longitudinal momentum fraction of the quark, one recovers the 
standard formula ^ 

dP = dzdlna. 

27r 1- z 

A further assumption that /? -C 1 leads to the double logarithmic approximation. 

Like the standard on-shell splitting function, eq. (|T5|) should be understood as an 
evolution kernel, so that multiple branching is described by iteration. This is achieved by 

DGE m]. 
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4 Gauge choice and interference 



Before proceeding with the apphcation of eq. (0), it is worth while to have a second look 
at the matrix element in order to appreciate the significance of the gauge choice. While 
the following discussion is quite general, to be concrete we will keep on using the example 
of e+e" — ^ qqg. Considering the hadronic tensor Hi^p, the rest of the amplitude in 
is J^r = l'^u{p), so the process-dependent trace in (0) becomes 



J2 Tr {^rMr^} = Tr {2/7 W} 



(16) 



If the angle with respect the electron beam is integrated over, as in the total fragmentation 
function case, it is sufficient to project H^p on —g'^''. Then one can replace Tr {]/7'^^7''} 
by 8pp. The light-cone gauge result for the squared matrix element, based on gluon 
attachment to the p quark only, is therefore 



light— cone 



ao 



-A 



1 + 2/3 + 2 



(17) 



(a + A)2 a + X (3 

where ao = SCpOs, but since the gluon attaches also to the antiquark p with an amplitude 



M 



^ {k+pf 



u^^\p)Y{^ + ^)Yu^'\p), 



;i8) 



the full squared matrix element is 



ao 



a/3 -A l + a + /3 + 2A 



I3a — A 



[a + Xf 



(a + A)(/3 + A) 



(/5 + A)2 
(19) 

This result is of course gauge- invariant. Let us now examine the origin of the various 
terms, in particular the singular terms for a + A — >■ 0, in different gauges. Let us compare 
the Feynman gauge and the light-cone gauge. In the Feynman gauge, M.M.'^ is the first 
term in (|19|) and the interference is responsible for middle term there. Thus in this 
gauge only part of the singularity is captured by a calculation {M.M.'^) that ignores the 
possibility of gluon attachment to the rest of the process. On the other hand in the light- 
cone gauge the entire a + A — > singularity is captured by (|l^ while the interference 
term is 



ao 



^ l + a + /3 + 2A ^ 1 + 13 



(20) 



(a + A)(/3 + A) ~{a + X)f3_ 

which is not singular. This shows that light-cone gauge is indeed useful in identifying the 
a -|- A — > singular terms. 

It is also interesting to note the different source of the singularity of the squared 
matrix element in the double logarithmic approximation. In the Feynman gauge the 
singularity originates in the interference term: 1/a appears because of the singularity of 
the propagator of the quark with momentum p + k, and l/j3 from the one with momentum 
p + k. In the axial gauge the singularity originates only from the diagram where the gluon 
interacts with p, so 1/a has the same source - the singularity of the propagator. On the 
other hand, 1//3 is just built into the gauge (H). 



9 



5 Dressing the gluon 



Let us demonstrate the application of the off-shell gluon splitting function for DGE. We 
begin with the simplest example: the case of the massles^ quark fragmentation in e~^e~ 
annihilation. The fragmentation functions were analysed in the dispersive approach by 
Dasgupta and Webber [32], who concentrated on the determination of the x dependence 
of the power corrections. Our results for SDG characteristic function in the limit x — > 1 
are consistent with their calculation. The main difference between our approach and 



that of ||2^, ^ is that in the latter the power corrections are additive to the coefficient 
function, while here (similar ideas were raised before, see |]18|-|22[, p6| , |27[] , and p3[]), 
as we shall see in the next section, the SDG calculation is considered as the kernel of 
exponentiation, so the power corrections will appear as a function that is convoluted with 
the resummed coefficient function. 

The fragmentation functions are defined by 



da 



a dxd cos 9 8 



-(1 + cos 



2 Frix, Q^) + ^ sin^ 9 Fl{x, Q^) + ^ cos ^ Fa{x, Q'), (21) 



where x = 2qph/Q'^- Here q is the electroweak gauge boson momentum (g^ = > 0), 
Ph is the hadron momentum and 9 is the angle of the hadron with respect to the e~ 
beam; Ft{x, Q^), Fl{x, Q^) and -^4(0;, Q^) correspond to the transverse, longitudinal and 
asymmetric fragmentation functions, respectively. Fa{x, Q^) vanishes for a pure electro- 
magnetic interaction. The standard factorization technique allows one to compute the 
coefficient function Cp{z,Q'^) by identifying the outgoing hadron as an outgoing par- 
ton i. The physical fragmentation function is then written as a convolution with the 
non-perturbative function Z)(x,Q^), 



Fp{x,Q')=J2 f'-CUz,Q')D,{x/z,Q'), 

i ^ 

where i can be a quark, an antiquark or a gluon. In moment space, 

Fp(x, Q^)x''-Hx = Fp{N, Q^), 

Jo 

the convolution translates into a product 

Fp{N,Q') = Y^CUN,Q')D,iN,Q'). 



(22) 



(23) 



We will now show that the log-enhanced terms in the quark fragmentation coefficient 
function C^(a;, Q^), which are associated with a single gluon emission, can be directly 
calculated from ([T5|) . In the next section this result will be used to calculate the multiple 
gluon (resummed) coefficient function through DGE. 

The gluon virtuality e has a crucial role in the calculation of the single gluon emission 
cross-section to all orders. For the calculation to be correct to any logarithmic accuracy 



perturbative analysis at NLL was recently completed for the case of heavy quark fragmentation. 
Here we discuss only the massless case. 
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(even in the simplest case, i.e. the large-Nj hmit), an integral over the running coupling 
must be performed with the correct argument - the gluon virtuality. This is the physical 
scale of the interaction. Using the dispersive approach in renormalon calculus, the integral 
over the running coupling is equivalent to the sum of the perturbative series in the large- 

limit. Such a renormalon calculation can be performed within the approximation 
considered in the derivation of (P^, where only the singular terms in the matrix element 
were kept. This approximation amounts to keeping all the logs in the differential cross- 
section, while neglecting terms that are explicitly suppressed by 1 — x. 

The total quark fragmentation function F{x, Q'^) = Ft{x, Q"^) + Fl{x, Q"^) is certainly 
sensitive to soft or collinear emission - we saw that the squared matrix element (|l^ (or 
(|19D) is singular for a + A — ^ 0. The quark that emits the gluon (in our gauge) is p. We 
therefore identify ph as p, getting 

x=^ = ^±^ (24) 

We see that the singularity in (^) is inversely proportional to 1 — a; ~ (a + A)/(l + /3). 
Using this variable, the differential gluon emission probability is 



dP CpOis 



dx 271 



dp 



1 ^ 1 p^ + 2p + 2 



x)2 (l+/3)2 1-x f3{l + (3y 



(25) 



where e = m?/Q'^. In the integration over (3, the limit that contributes to log-enhanced 
terms is the limit where the gluon transverse momentum 7 vanishes. This translates into 

f3 = \/a = e/{l-x-e). (26) 

We therefore integrate over [3 and substitute (p6D, getting 

dP Cpas 



dx 2tt 



H^.e) (27) 



witl 



M - , /'l-a^X 3 1 e 1 
H^.^)\lo. = In - + n ^ + o 7^ ^- (2^ 



1-x V e / 2 1-a; (1 - x)^ 2(l-x)3' 
The physical region is e < 1 — x, where a and (3 are positive; J-'{x, e) is the off-shell gluon 



characteristic function [^9], |3^, ^ for the quark fragmentation function. The subscript 
'log' in ( P5| ) is meant to remind us of the approximation: by taking only the 1 — x — > 
singular terms in the matrix element, we control only log-enhanced terms in the pertur- 
bative coefficients. The full characteristic function (see e.g. [^) contains some additional 
terms, which generate non-logarithmic terms in the perturbative series. 



^Here we dropped sub-leading terms (e.g. 0(e^/(l — a^)^)), which go beyond our approximation and 
contribute only to non-logarithmic terms in the perturbative coefficients. 
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In the dispersive approach |^ (see also [2B|, |T^, 0) the single off-shell gluon result 
is promoted to an infinite perturbative-sum associated with a single dressed gluon by 
replacing the coupling with a dispersive integral 



2(3o Jo 



1-^' de ■ 



2/3o Jo 



1-^- de 



(29) 



where A{k'^) = fiQasik"^) /t^ ■ The two integrals in ( pQ]) are related by integration by parts: 
jF(x, e) = —edJ-'{x, e)/de , and the function p(/i^) is identified as the time- like discontinuity 
of the coupling, 



p(/i2 



— Disc|i(-u2)| = — 
27ri I V ;/ 27rz 



A -p" + a -A -p"-zA 



(30) 



The "time-like coupling" A^g{fi^) in (^91) obeys dA^^i.!^^) / djJi^ = p{fi^)- For example, in 
the one-loop case A{k'^) = l/(lnA;^/A^) and 

^cff (/U ) = 77 arctan - In . 

2 vr yvr A^y 

We emphasize that (^) is considered just as a perturbative sum, and not as a non- 
perturbative definition^. 

The bar over the coupling indicates a specific renormalization scheme. In the large-/3o 
hmit (large-iV/ limit) A{k^) is related to the MS coupling Aj^{k'^) = Poof^ {k'^) / n by 
A{k'^) = A-^{k'^)/ (^1 — I^Mg (A;^)^. In the general case, going beyond this limit requires 
a skeleton expansion |T^, generalizing the Abelian one. However, in the context of our 
approximation, where only the singularity of the splitting function matters, a skeleton ex- 
pansion is not required. The appropriate coupling is the "gluon bremsstrahlung" effective 
charge PO, 123] . Fixing A by 



Ams {k'^ 



1 



+ 



(31) 



the splitting function (^) is correct to next-to-leading order, as far as the singular terms 
are concerned. The origin of the "gluon bremsstrahlung" effective charge can be under- 
stood by the formulation of the Sudakov resummation through an evolution equation for 
Wilson-line operators. In this language this effective charge appears as the cusp anoma- 
lous dimension ^ . 

In the quark fragmentation case (^81), we have 



j^{x, e] 



log 



1 — X 



[l-xf 



(32) 



Note that J^{x, e = 0) 7^ so the integral (^) diverges for vanishing gluon virtuality. This 
is a direct consequence of the fact that the fragmentation function is not coUinear-safe. 

^Such a definition differs from ttie Borel-sum by power-corrections wtiicti are not associated witti tlie 
infrared dynamics p^ . 
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Factorization is required. We shall make the following identification: Q"^) = 

6{1 — x) + C|nQ(x, Q^) + . . . where the eUipses stand for corrections that are sub-leading 
by powers of either Pq or 1 — x and 



_Cf \ 2 



-S{x) A.,(eQ2) + I - {t{x, e) - S{x)) A^^ieQ') 



1 — X Ji 



1 e 
1-^ de 



e Jo 




(33) 



+ 



XI 



Here S{x) = 2/(1 — x) is the coefficient of ln(l/e) in J-'{x,e). The subtracted term is 
absorbed into the non-perturbative function Dg{x/ z,Q'^). This means that its evolution 
is controlled by S{x), i.e. by Ci?(Q;s/7r)/(l — x). With the definition (|33|) the coefficient 
function has a well-defined perturbative expansion. Expanding the coupling A^g^eQ'^) in 
a fixed-scale coupling one can write the corresponding perturbative coefficients in terms 
of the log-moments: 



hJx) = 



1-^- de / 1 
— m - 
1 eve 



x] 



n + 1 



-L 



Six) 



n+l 



1-^- dt /, 1 
— m - 
e V e 



Six 



(34) 



e'^r {n + 1, L) + 2~"~^ e'^^F (n + 1, 2L) 



where L = lnl/(l — x). Since e^r{n + 1,L) = Y^^^q n\ / j\ , the expression in the square 
brackets is just a polynomial of order n + l in L. Note, however, that since the sub-leading 
logs increase factorially, the sum is not well approximated by the leading log term even 



at rather large L (see P3|). 

Contrary to the total fragmentation function, the longitudinal fragmentation function 
F^iXjQ"^) is insensitive to soft or collinear emission (at the logarithmic level): projecting 



(0) onto elel 



where 




IpI V Q 

is the polarization along the direction of p, the contribution vanishes. It follows that 
in the electromagnetic case the log-enhanced terms in the total fragmentation function 
all originate in the transverse fragmentation function (see [p2 , 



Exactly the same 



log-enhanced terms appear in the parity-violating asymmetric fragmentation function. 



6 Dressed gluon exponentiation 

We saw that factorization of the squared matrix element, which is the basis for exponen- 
tiation, holds also for an off-shell gluon. Exponentiation can be derived, as in Q , from 
an evolution equation, or, as in , directly from the probabilistic interpretation of the 
SDG cross-section based on the approximation of independent emission. 

The factorization of the squared matrix element suggests that the entire log-enhanced 
SDG sum exponentiates. In reality straightforward exponentiation holds only up to a 
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certain logarithmic accuracy, because gluons are correlated through the phase-space in- 
tegration. Momentum conservation translate into a factorized form in Laplace space (see 
e.g. |2^, ^). To some logarithmic accuracy one can then perform the integration over the 
gluon momenta as if other gluons had not been emitted. At the same time exponentiation 
of the entire SDG sum can be regarded as the leading contribution to the exponent in the 
large- /?o limit. In this sense the expansion in powers of l//5o in the exponent replaces the 
standard expansion in powers of the log. 

Being aware of this we now proceed to write down the explicit exponentiation formula 
for the fragmentation function case. We work in the approximation where the invariant 
mass of the jet, Q^{1 — x), is additive with respect to multiple gluon emission. With this 
assumption we ignore correlations between the gluons, which would modify the exponen- 
tiation beyond NLL accuracy. The resummed expression for C"'(x, Q^) is then given by 
the inverse Mellin transform of Cq{N,Q'^)\^^^, where 



C,{N,Q' 



exp 



dx \x 



1) 



(35) 



Here we drop an x-independent factor associated with the virtual corrections to the hard 
cross-section (see e.g. [0). 

A convenient way to deal with perturbation theory to all orders is Borel summation. 



Starting with the scheme-invariant Borel representation of the coupling |p6||-p8 

A{e) = du exp (-U In ^ / K^) Ab{u) 



and taking the time-like discontinuity we get 



du exp (— wlnQ^/A' 



smvr-u -r 



TTU 



Ab{u). 



(36) 



(37) 



Using this Borel representation of A^gi^eQ"^), and integrating (^) over e, one obtains a 
Borel integral representing the SDG perturbative sum: 

poo 

C|dg(^'<5^) = 7T^ / duBsDG{u,x) exp(-u\nQ'^/A'^ 
with the following Borel function 



smTTU - , , 
Ab{u) 

TTU 



{31 



-1 
1 — X 



2 ((1 - x)-" - 1) ^ (1 - xy ^ (1 - a;)~" 



u 



1-u 



2-u 



(39) 



Next, the exponent in Mellin space In Cq{N, Q^) |^^^ is obtained from ( |35D integrating over x 



\nCqiN,Q^) = ^ r duB^{u) exp (-m In Q VA^) ^-^^ Ab{u) 
2po Jo ^ ^ nu 



(40) 



with0 



Bn(u) = — 



u 



e«i'^^r(-iz) + - + 7£; + lnAr 
u 



1 



+ 



u 



u 



-u) + - 
U. 



^Upon integration V{—u) appears with a factor r(iV)/r(iV 
(A'' ^ 1) this factor can be replaced by e"'"^. 



(41) 

In the approximation considered 
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The Borel function (^T]) contains valuable information on both the perturbative and 
the non-perturbative levels. 

Sudakov logs: To get an explicit perturbative expansion of the exponent in Mellin 



space one can perform the Borel integral term by term, similarly to what was done in p3 



for the case of the thrust distribution. For simplicity we work with a one-loop running 
coupling y4(Q^). It is straightforward to rewrite the sum as 

/~i oo 

lnC,(iV,Q2) =_^Y.Am''' A(A(Q2)lniV) , (42) 
where the functions fk{C) sum all powers of ^ = A{Q'^) \n.N . They are given by 

(fc-l)/2 ( 2\l 

M0= E -7^^rTiv9!\0MUk), (43) 

where gl!\^) is the following hypergeometric function 

9!\0 = r ([1, A; - 2 + no], [1 + no], (44) 



with no = max(0, 21 + 3 — k) and 



fc-2 



M,oi{k) = 2 Ck-2 + E (1 + 2"""') Cfc-3-m, (45) 



m=0 



where the numbers are defined by T{—z) = — c^z''. This means that c_i = 1, 

Co = 1e, Ci = 7r^/12 + 7|;/2, etc. Higher contain higher Q numbers {i < k + 1), yet 
numerically Ck {k > 1) are all close to 1. For k = 1,2 the functions are 

= 2 (e + (l-01n(l-0) (46) 
/2(0 = -(l + 2^E)\na-0 



while for k > 3 they are characterized by increasing coefficients as well as increasing order 
pole singularities at ^ = 1, 

/3(0 = 0.804/(1 -0 + 3.29 e 

/4(0 = 0.779/(1 -0' 

/5(0 = 2.32/(1 -0'-9.74e 

/6(0 = 6.12/(1-0^ .47. 

/7(0 = 23.8/(1-0' + 45.78 e ^ ' 

/8(0 = 120.9/(1-0' 

/9(0 = 721.54/(1 -0^- 263.57^ 

/io(0 = 5043.48/(1 -0^ 

We see that sub-leading logs are characterized by an explicit factorial growth as well as 
by an increasing singularity at the end of the perturbative region: 

e^^ln-^ = emax, (48) 
TT 1 — X 
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where ^max = 1. A similar structure was identified in the case of the thrust (or jet mass) 



distribution in [g3|. The two cases are compared in the next section. As discussed in p3[ 
the enhancement of sub-leading logs implies that an approximation based on a fixed 
logarithmic accuracy in the exponent has a small range of validity, namely C, ^ ^max- 
The resummed exponent (^OD is valid much closer to the strict limit of applicability of 
perturbation theory ^ ^max), provided appropriate power corrections are included. 

Power corrections: The form of the power corrections can be deduced from the 
renormalon ambiguity. They are expected to be additive at the level of (^Op and therefore 
they exponentiate together with the perturbative sum, so that 

C,{N, Q') = C,{N, Q') C^^iN, Q'). (49) 

res 

This implies that the correction can be written as a convolution in x space. The Borel 
singularity of the fragmentation function exponent in the large-/5o limit can be read off 
eq. d^), taking into account the attenuation by the sin(7r'u) factor in (^0]). The result 
turns out to be very simple: a pole at z = 1 with a residue — and second pole at z = 2 
with a residue —N'^/4. The non-perturbative correction C^^(A^, Q^) is therefore expected 
to be 

where lOn are dimensionless constants. Since Cq{N, Q"^) multiplies a non-perturbative 
distribution Dg{N,Q'^), eq. ( ^0]) leads to a meaningful prediction only if it dominates 
the behaviour in the threshold region. This possibility should be tested experimentally. 
As discussed in the next section, the situation is different for collinear- and infrared-safe 
quantities where no non-perturbative distribution is needed a priori. 



7 Jet mass in 



e+e 



annihilation 



One of the classical examples where soft and collinear gluon radiation is important is the 
case of event-shape variables in e~^e~ annihilation. As opposed to the fragmentation func- 
tion case, these observables are infrared- and coUinear-safe, so the perturbative calculation 
does not require factorization. The enhanced sensitivity of these observables to large angle 



soft emission is reflected both perturbatively - as large logs and non-perturbatively 
- as power corrections [|T3|, ^ (see also |29|, |3|]-[||] and lT|]-[g). 

Here we consider event-shape variables that are related to the jet mass distribution, 
such as the thrust and the heavy jet mass. The thrust distribution was analysed in detail, 
using DGE, in Analysis of the heavy jet mass by the same method will appear 

soon [Q. The starting point in ||2^ was the full SDG characteristic function We 
shall now see how the relevant part of the characteristic function, the one that generates 
the logs, can be computed directly from the generalized off-shell splitting function. This 
will be followed by a brief description of the results of for comparison with other cases 
studied here. 

Consider a single off-shell gluon emission, which is either soft or collinear (two-jet 
kinematics). Assuming that the gluon is in the hemisphere of the quark p, it is p that 
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sets the thrust axis. The p jet mass is zero and the p jet mass, which is also 1 — thrust, 
is simply 



P 



'"^H (P + 2pA; + A;^ (a + A) 2pp 



a + X 



l + a + P + \ 



1 — X. 



(51) 



Thus calculation of the SDG characteristic function (neglecting non-inclusive effects ^9 
pT^ , |43| , |23[| , associated with gluon decay into opposite hemispheres) is just identical to 



that of the total fragmentation function. The result is 
by p: 



with the replacement of 1 



■X 



Hp. 



(52) 



- 2 

log p p^ 

This result is confirmed by the explicit calculation of the full characteristic function 



followed by the identification of terms that contribute to logs |23 



A major difference between jet mass distribution and the fragmentation function case 
is in the large- angle phase-space limit: we assumed that the gluon is in the hemisphere of 
the quark p. This is equivalent to (3 > a. Along the zero transverse momentum boundary 
of phase space a(3 = A, this condition translates into (3 > -\/A. It follows that the lower 
integration limit over e is e = p^, which guarantees infrared- safety. The upper integration 
limit is e = p, as for the fragmentation function case. The SDG differential cross-section 
is therefore 



a dp 



-f / A,«.(eg2)^(p,e), 



(53) 



where J-'(p, e) is given by (m) . 

Proceeding along the lines of the previous section we first use (^7|) to write a Borel 
representation of the SDG cross-section. The Borel function is p3[| 



Bsog{u,P) 



2 / , r 

— exp 2u In - 
« V P. 



Next, the exponent in Laplace space p3, p3 



2 

u 

is 



1-u 



Jo a dp ^^^V 



The corresponding Borel function is given by 



B,{u) 



dp 



p 



— e 

u 



2u Ini 



1 1 
+ 



U 



2u In u 



r(-2n) + 



u 



2u 



1-u 



1-u 
1 



u 



^ exp \ u In - 

2-u) p, 



-vp 



u In - 



U 



1) dp. 



(e-- - 1) 



(54) 



(55) 



(56) 



u In V 



Ti-u) + 



u 



where non-logarithmic terms were neglected. It is interesting to compare this result with 



the fragmentation function case (^Of ) and (|41|) , where only the coUinear limit e = 1 — x = p 
plays a role. The terms asociated with the coUinear limit are common. Equation ( p6D 
contains in addition a term that corresponds to large-angle emission e = p^. The latter 
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has a double argument u — ^ 2-u, indicating larger perturbative coefficients as well as 



enhanced power-corrections . 

Equation (^) summarizes the DGE result for the jet-mass distribution. Performing 
the Borel integral as in (EOl) and the inverse Laplace transform, this result can be directly 



used for phenomenological analysis of the thrust and heavy jet mass distribution |23, 44 



A detailed analysis of the perturbative structure of the exponent (sub-leading logs) and of 
the expected power corrections was performed in 123] . The main results are the following. 



Sudakov logs: Performing the Borel integral term by term with the one-loop running 



coupling we can write the exponent In Jj,((5^) as in eq. (^21) with 



(fc-l)/2 / o.i 



fkiO= E 



(57) 



7-0 (2^ + 1^ 

where ^ = A{Q'^) Inz/ and gk\0 is given by (^). Mcoi(fc) (given by (^)) and 

M^k) = -2''-'ck-2 (58) 

correspond to the collinear and large-angle contributions to the Borel function, respec- 
tively. Explicitly, this yields 



MO 



2(1-0 ln(l-0-(l 
-2 7E (ln(l-0-ln(l 



2 ln(l-20 

2 0)- ?ln(l-0 



(59) 



and 



/3(0 = 

m) - 

/5(0 = 

m) - 

friO - 
fsiO - 

f9{o - 

/io(0 



0.804/(1 -0 -2.32/(1 -20 
0.779/(1-0^-2.68/(1 -20^ 
2.32/(1-0^-5.00/(1 -20^ 
6.12/(1-0^ - 15.32/(1 -20^ 
23.8/(1 - 0^ - 61.12/(1 - 20^ 
120.9/(1 - 0^ - 305.52/(1 - 20^ 
721.54/(1 - 0^ - 1833.55/(1 - 20^ 
= 5043.48/(1 - 0® - 12834/(1 - 20®, 



The factorially increasing coefficients of /fc(0 ^^id the enhanced singularity at ^ = 1/2 
imply that a fixed logarithmic accuracy calculation holds only in a very restricted region 
^ <^ 1/2. Otherwise, these sub-leading logs must be resummed |2^. Having performed 
such resummation and having included the appropriate power corrections (see below) the 
region of applicability of the result can be extended up to ^ < ^max = 1/2- We saw that 
in the absence of large angle soft emission sensitivity (the fragmentation function case) 
this region is larger: ^ < ^^ax = 1- 

Power corrections: The ambiguity of the Borel integral indicates specific power 
corrections Taking into account the singularities of ( ^6|) and the factor simru/iru of 
eq. (|37|), the Borel integrand (in the large-/3o limit) is singular at half integers m = |, |, | . . ., 
and at M = 1,2. The ffist type of singularity is a result of the large angle soft emission 
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(the first term in (^)). It translates to power corrections of tlie form ~ [yjQY" where n 
is an odd integer, 

2i-l' 
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The large angle power corrections are important when p 
can be resummed into a non-perturbative shape function 
inverse Laplace transform of (pO). Thus, the calculation of the exponent in the large-/3o 



Xjv approaches KjQ. They 
23|| , corresponding to the 



limit by DGE yields highly non-trivial information on the shape function: its even central 
moments (controlling corrections ~ iylQY' with even n) are suppressed |33|. The second 
type of singularity is a result of coUinear emission. As in the fragmentation function 
case, it translates into power corrections of the form ~ l/{Q'^p) and ~ l/{Q^p^). These 
corrections can be neglected in the region pQ/ ^ A. 



8 Deep inelastic structure functions 

Deep inelastic structure functions differ from the previous examples by the fact that they 
have a systematic expansion in the powers of the hard scale - the twist expansion [^, 
In the region of interest, Xbj = Q^/(2Pg) — >■ 1 (here P is the nucleon momentum and q 
is the virtual photon momentum (— = Q^)) the twist expansion breaks down and power 
corrections of the form 1/Q^"(1 — xbj)" become dominant. In this respect, the situation 
is analogous to the fragmentation function case. 

The renormalon approach was used to study the xbj dependence of power correc- 
tions (see ^"^^ HD)- The correspondence between the renormalon approach and 



the twist expansion was addressed in |T^ in the specific case of the longitudinal structure 
function F^. It was shown there that the infrared renormalon ambiguity of the twist-two 
coefficient functions cancels against the ambiguity of the twist-four matrix elements, am- 
biguity that results from the ultraviolet quadratic divergence of these matrix elements. 
The case of the transverse structure function F2 is more interesting than that of F^, owing 



to the threshold region I^Sj. The matching of the renormalon ambiguity and the twist 
expansion in this case is discussed in [p^]. 

Contrary to the e^e~ examples analysed above, in the case of deep inelastic scattering, 
both incoming and outgoing quarks can be the source of gluon bremsstrahlung. Never- 
theless, in the approximation considered (and in the appropriate gauge) one can describe 
the radiation as if it were associated with the final-state quark alone. Indeed, we will see 
that the SDG characteristic function, which generates the log-enhanced part of the F2 
coefficient function, is identical to that of the quark fragmentation function. 

Identifying the incoming quark momentum P with p ( "— " direction) and fixing the 
axial gauge ^4+ = we consider gluon emission {k) off the outgoing quark, which carries 
momentum p in the direction. In this gauge the gluon coupling to the incoming quark 
does not give rise to singular terms. Momentum conservation implies that p + q = p -\-k. 
The standard xbj, expressed in terms of the light-cone parameters (|^) is 

XBj = QV(2Pg) = QV(2w) = !-(« + A)/(l + /3). (61) 
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The matrix element for emission off the outgoing quark p is (|^). The phase-space 
integral is 



/ ^ " "^'^ ^ (^') + (62) 

I f dad/3 ^. ..^ 
^8^J + ^ (" - (1 - ^Bj - ea;Bj)(l + P)) , 

where e = k'^/Q'^ = A/(a;Bj(l + P))- Expressing in terms of xbj and /?, we find that 
the partonic differential cross-section for a single gluon emission coincides with (^) upon 
replacing 1 — a; by 1 — xbj- As in the e^e~ case, the integration over j3 is bounded from 
below by the positivity of the squared transverse momentum. For vanishing transverse 
momentum a(3 = X one finds 

P=^^ (63) 

1 — XBj — exBj 1 — xbj — e 



Integrating the gluon emission probability over /? and substituting the limit (|5BD, one gets 
the following characteristic function 



2 e 

^ix,e) = - r^, (64) 

log 1 - XBj (1 - XBjY (1 - XBj)^ 

which is identical to the fragmentation function case (0). This result can be compared 
with the full characteristic function computed in |2^. The additional terms in the 
latter do not contribute to log-enhanced terms in the perturbative coefficients. In the 
approximation considered, the integration over the gluon virtuality has the same range 
as in the fragmentation function case, i.e. e < 1 — xbj- Also here there is no restriction 
from below (the coefficient functions are not collinear-safe) so that factorization into a 
non-perturbative parton distribution must be used. It follows that the structure of the 
exponent in the large-/5o limit is just identical to that of the fragmentation function case, 
as summarized by eqs. and (p]). As discussed at the end of section 6, this structure 
has definite implications in both the perturbativef] and non-perturbative levels. The fact 
that the exponent is identical to the fragmentation function case suggests that it is the 
general properties of the quark jet evolution that control the power corrections and not 
the particular process from which the jet emerges. 



9 Initial-state radiation in the Drell— Yan process 

The classical example where initial-state radiation determines the cross-section near thresh- 
old is the case of the Drell- Yan process, the inclusive cross-section for lepton pair pro- 
duction in hard hadronic collisions. Resummation of the leading and next-to-leading logs 



"The resummation of deep inelastic coefficient function was recently performed |53 to NNLL accuracy. 
Our results are consistent with those calculations. 
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for this process in the hmit where x = /{p + pY approaches 1 was 
and the associated power corrections were discussed in p^, E6|, p] 



performed ^, [50 



The factorization property of soft and colhnear gluons holds for initial-state radiation 
exactly as for final-state radiation. In particular, at the level of the squared matrix 
element, the emission of soft or colhnear gluons appears as an overall factor, multiplying 
the cross-section for no emission, as in eq. (^. However, as we shall see below, the way 
the singularity of the propagator is related to the physical parameter that controls the 
phase space, 1 — x, is fundamentally different and consequently so is the structure of the 
exponent. Although the leading and next-to-leading logs are similar to those of the jet 
mass (or the thrust) distribution in e^e~ annihilation, the characteristic function and 
thus the sub-leading logs and the power corrections 



26 are different. 



Consider now the case of quark (p) antiquark [p] annihilation into an off-shell boson q 
(a photon in Drell-Yan), where the quarks originate in the colliding hadrons. In the 
perturbative calculation, we begin with on-shell partons p"^ = p^ = 0. Choosing the frame 
such that p is in the "+" direction and the gauge = 0, so that the radiation decouples 
from p, we start by calculating the squared matrix element for a single gluon emission 
off the quark p. The gluon momentum k is decomposed as in (^. Using (|TD|) and the 
following phase-space integral 



/ 7 — -:r 5 ik"^ — d'^qd (q^ — 2ppx] 5^{q + k — p — p) = [ 5(1 — x — [5 — a + X), 

J (27r)'^ V / V / J gy;-^ 



dadf3 



(2vr) 

we find that the partonic differential cross-section for a single gluon emission is 



1 da Cpa, 
a dx 



2tt Jpi 



P2 



d(3 



-A 



{i-x-py- 



(1-/9) + 



/32 - 2/3 + 2 



l-x- (3 



/9 



(65) 



The integration over (3 should now be performed for fixed im? and s = 2pp, thus for 
a fixed A. As in the previous cases, the relevant integration limit is deduced from the 
condition that the gluon transverse momentum vanishes, i.e. af3 = A. Substituting the 
relation a = l — x — (3 + X, one obtains a quadratic equation for f3 whose solutions are 
A ± A)/2, where A = J{1 - xf - 2A(1 + x) + A^. Integrating {E^) one 



/3i 2 = (1 - a; + A ± A)/2, where A = ^(1 
obtains the following characteristic function: 



^(x,A)| 



log 



In 



X 



[l-x-X + A){l-x + X + A) 
;i-x-A- A)(l-x + A- A)' 



(66) 



where, as in (|28|), a factor of CfO!s/{2tt) was extracted. Since in this case only the small 
f3 region contributes to log-enhanced terms, the only relevant term in (|65|) is 



1 da 
a dx 



2tx 



di3 



X 



(67) 



This simple expression is sufficient for the calculation of (|66D. 
Taking the logarithmic derivative with respect to A yields 



jF(x, A) 



log 



^(1 -x)2 -2A(l + x) + A2' 



(6^ 
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Comparing this characteristic function to the case of e^e~ annihilation or (^) and 
deep-inelastic scattering (|6^, it is clear that the structure of the perturbative expansion 
and the power corrections will be completely different. This should be contrasted with 
the picture one might get based on the NLL resummation formula: in the so-called deep 
inelastic factorization scheme 
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the latter is identical to that of the thrust. In fact, 
the similarity of the leading logs is superficial. There is a fundamental difference between 
the first set of examples, where the radiation originates in an outgoing quark (an evolving 
quark jet) with a constrained invariant mass, and the Drell-Yan case, where the radiation 
originates in the mcominf/ quarks and the total energy {a+(3) in the final state is restricted. 
As we shall see this difference is encoded into the structure of the exponent making an 
impact on both the sub-leading Sudakov logs and the power corrections. 

We would like to express the SDG cross-section in the form of a Borel sum, as in (^). 
We should therefore calculate / (iA A^"^^jF(x, A). As in the cases of the fragmentation 
functions and deep inelastic structure functions, the lower limit is zero so the integral 
diverges. Factorization must be used to obtain a well-defined Borel function. Using the 
same subtraction prescription as in (^3]) we get 



(l-v^) 



1 — X Jl 



ciAA-"-i 



(69) 



(l-x)2-2A(l+x) + A2 



The upper integration limit was deduced from the condition /?2 < yielding A(Ama 
and therefore Amax = (1 — Since the integration over A is restricted to A ^ 1 



the logarithmic approximation one can replace the integration limit by A 
and approximate the A-dependent denominator in 



1 



X] 



Iby 



X 



2/4 



2 - 4A. We thus find 



Bsuc{x,u) 



1 



\ — X u 



1 - 



rfi -u) /i-x 



r(i-«) 



-2u 



(70) 



Finally, the Borel function corresponding to the exponent in Mellin space is 

Bn{u) = £dx (x^-i - l) B,^o{x,u) = 2 (e'"^"^ - l) T{-uy - ^IniV- (71) 

As far as the leading terms at large are concerned, this result agrees with that of 
ref. 1 2^, which used the full matrix element. 

Sudakov logs: Equation (^) allows us to examine the structure of the exponent to 
a fixed logarithmic accuracy, as in (^). Here fkiO are given by 



(fc-l)/2 



fk{0= E 



{21 



(72) 



where gk\C,) is defined in ( ^4|) and 



k-2l-2 

Min(/c,/)=2 Ck-2l-3-jCj. 

i=-i 



(73) 



22 



Explicitly, this yields 



/i(0 = 2(1-20 ln(l- 20 +4e 

MO = -47ln(l-20 (74) 

and 

/3(0 = 1.33/(1- 20 + 6.58e 

/4(0 = 2.12/(1 -20' 

/5(0 = 4.00/(1 -2 0=^- 19.48 e 

/6(0 = 11.59/(1 -20' 

/7(0 = 48.42/(1 -20' + 91.56 e 

/8(0 = 238.80/(1 -2 0' 

/9(0 = 1438.66/(1 - 20' - 527.14^ 

/io(0 = 10078.0/(1- 2 0^ 

As in the previous examples, sub-leading logs are enhanced by factorially increasing co- 
efficients as well as an increasing singularity at ^ = ^max = 1/2. 

Power corrections: The power corrections in the x — > 1 region can be deduced 
from the singularity of the Borel transform of the exponent (|7TD . As in the previous 
examples, they are expected to exponentiate together with the perturbative sum. The 
singularities of the Borel integrand are located at integer values of u. Equation ( [Til) has 
double poles at all integers, but the sin(7rM)/7rM factor of (pT]) leaves only simple poles. 
As observed in [^, this singularity structure implies that the leading power correction 
at large (and not too large x) is 1/(5^(1 — x)^. Closer to x = 1 sub-leading power 
corrections of the form 1/{Q'^{1 — x)^)", where n is an integer, become important. The 
non-perturbative correction factor (in Mellin space) is 



CMN, Q') = exp 5: /\ ' u;„ -— . (75) 




10 Conclusions 

The factorization of soft and collinear gluon bremsstrahlung is one of the fundamental 
tools of perturbative QCD [jl|. In particular, it sets the basis for resummation of loga- 
rithms in many applications. In this paper we demonstrated that factorization holds also 
in the case of an off-shell gluon. The main advantage in keeping the gluon off-shell is that 
the virtuality provides the right physical scale for the coupling. This opens the possibility 
to resum a full set of radiative corrections, which are associated with dressing the emit- 
ted gluon (a renormalon sum). Keeping the gluon off-shell, we identified the appropriate 
kinematic approximation needed for the calculation of log-enhanced terms originating in 
either the collinear or the soft regions of phase space. The factorization formula led to 
the definition of a generalized splitting function that describes the emission probability of 
an off-shell gluon off a quark. As always, factorization results in exponentiation. Thus, 
DGE promotes the leading order calculation to an all-order sum in two respects: the first 
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by dressing the gluon and the second by taking into account multiple emission through 
exponentiation. 

In the class of problems analysed here, where a kinematic threshold imposes a stringent 
restriction on the real emission phase space, resummation is necessary even for a qualitative 
description of the differential cross-section. Currently available and future experimental 
data provide a challenge to QCD: quantitative predictions for differential cross-section are 
required. The threshold region turns out to be phenomenologically important in many 
cases: for example in e~^e~ annihilation the cross-section is large only in the two-jet region. 

The analysis of the exponent in the different examples considered (see also |2^) shows 
that sub-leading Sudakov logs are enhanced factorially with respect to the leading logs. 
This is a direct consequence of the integration over the running coupling, which reveals the 
presence of infrared renormalons. Thus because of the running coupling, a quantitative 
description of the threshold region, e.g. in the fragmentation function case Q'^{1 — x) ^ A 
or {Po<^s/ vr) L ^ 1 (where L = In 1/(1 — x)), requires the calculation of the exponent to any 
logarithmic accuracy. Resummation with a fixed logarithmic accuracy, e.g. NLL, applies 
only for (/?o«s/7r) L <C 1. Moreover, in the threshold region, non-perturbative corrections 
are particularly significant. As opposed to single-scale observables, non-perturbative cor- 
rections are not dominated by a leading power, but rather appear as a convolution of 
the perturbative distribution with some shape function of the relevant scale. 

We saw that the exponent calculated through DGE contains both perturbative and non- 
perturbative information. As in the case of single-scale observables JlSl, these two aspects 
cannot be separated and must be treated together by resummation and parametrization 
of power corrections. 

As any renormalon-based approach, DGE does not pretend to supply a field theoretic 
definition of the shape function or to determine the magnitude of the non-perturbative 
corrections. However, it strongly suggests: (1) the factorization of the non-perturbative 
corrections; (2) the scale on which the non-perturbative function depends, e.g. (1 — x)Q'^ 
in the fragmentation and deep inelastic cases and Q{1 — T) and Q{1 — x) in the thrust 
and Drell-Yan cases, respectively; and (3) the specific powers of this scale appearing in 
the exponent (in moment space), for example: odd powers of u/Q in the thrust, {N/Q'^)"', 
where n = 1 and 2, in both the fragmentation and deep inelastic cases, and even powers 
of N/Q in the Drell-Yan case. The first two properties can be deduced also from other 
considerations. The third is harder to access non-perturbatively. Here one really relies on 
the assumption that the dominant non-perturbative effects are the ones detected by the 
perturbative tools. 
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